Optimal quantum estimation in spin systems at criticality 
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It is a general fact that the coupling constant of an interacting many-body Hamiltonian do not 
correspond to any observable and one has to infer its value by an indirect measurement. For 
this purpose, quantum systems at criticality can be considered as a resource to improve the ul- 
timate quantum limits to precision of the estimation procedure. In this paper, we consider the 
one-dimensional quantum Ising model as a paradigmatic example of many-body system exhibiting 
criticality, and derive the optimal quantum estimator of the coupling constant varying size and 
temperature. We find the optimal external field, which maximizes the quantum Fisher information 
of the coupling constant, both for few spins and in the thermodynamic limit, and show that at the 
critical point a precision improvement of order L is achieved. We also show that the measurement 
of the total magnetization provides optimal estimation for couplings larger than a threshold value, 
which itself decreases with temperature. 

PACS numbers: 03.65.Ta, 05.70.Jk, 06.20.-f 



I. INTRODUCTION 

Acquiring information about a physical system in- 
volves observations and measurements, whose results are 
subjected to fluctuations, and one would like to elimi- 
nate or at least to minimize the corresponding errors. 
However, the precision of any measurement procedure is 
bounded by fundamental law of statistics and quantum 
mechanics, and in order to optimally estimate the value 
of some parameter, one has to exploit the tools provided 
by quantum estimation theory (QET) [l|. 

As a matter of fact, many quantities of interest do not 
correspond to quantum observables. Relevant examples 
are given by the entanglement or the purity of a quan- 
tum state [2| or the coupling constant of an interacting 
Hamiltonian. In these situations one needs to infer the 
value of the parameter through indirect measurements. 
For many-body quantum systems, changing the coupling 
constant drives the system into different phases and, in 
turn, this may be used to estimate the coupling itself. 
In particular, close to critical points, quantum states 
belonging to different phases should be distinguished 
more effectively than states belonging to the same phase 
[S, 4, 5, 6, 7, 8, 9]. Distinguishability is usually quanti- 
fied by fidelity between quantum states, i.e. overlap be- 
tween ground state wave functions. In turn, the fidelity 
approach to quantum phase transitions (QPT) has re- 
cently attracted much attention [3, [ll| since, differently 
from bipartite entanglement measure approach [l2l |. it 
considers the system as a whole, without resorting to bi- 
partitions. In estimating the value of a parameter, one 
is led to define the Fisher information which represents 
an infinitesimal distance among probability distributions, 
and gives the ultimate precision attainable by an estima- 
tor via the Cramer-Rao theorem. Its quantum counter- 
part, the quantum Fisher information (QFI), is related 
to the degree of statistical distinguishability of a quan- 



tum state from its neighbors and, in fact, it turns out to 
be prop ortional to Bures metric between quantum states 

[HE [H, [13, [H, [H, [13, [H. 

As noticed in [ij] one can exploit the geometrical the- 
ory of quantum estimation to derive the ultimate quan- 
tum bounds to the precision of any estimation proce- 
dure, and the fidelity approach to QPTs to find working 
regimes achieving those bounds. Indeed, precision may 
be largely enhanced at the critical points in comparison 
to the regular ones. Here we show that the general idea 
advocated in |13| can be successfully implemented in sys- 
tems of interest for quantum information processing. To 
this aim we address a paradigmatic example of many- 
body system exhibiting a (zero temperature) QPT: the 
one-dimensional Ising model with a transverse magnetic 
field. 

In most physical situations, some parameters of the 
Hamiltonian, e.g. the coupling constant, are unaccessi- 
ble, whereas others may be tuned with reasonable control 
by the experimenter {e.g. external field). Therefore, the 
idea is to tune the controllable parameters in order to 
maximize the QFI and thus the distinguishability and 
the estimation precision. In doing this we consider the 
system both at zero and finite temperature, and fully ex- 
ploit QET to derive the optimal quantum measurement 
for the unobservable coupling constant in terms of the 
symmetric logarithmic derivative. In the thermodynamic 
limit we find that optimal estimation is achieved tuning 
the field at the critical value, in accordance with [l3j . 
whereas at finite size L, the request of maximum QFI 
defines a pseudo-critical point which scales to the proper 
critical point as L goes to infinity. In turn, a precision 
improvement of order L may be achieved with respect to 
the non critical case. 

The optimal measurement arising from the present 
QET approach may be not achievable with current tech- 
nology. Therefore, having in mind a practical implemen- 



tation, we consider estimators based on feasible detec- 
tion schemes, and show, for systems of few spins, that 
the measurement of the total magnetization allows for 
estimation of the coupling constant with precision at the 
ultimate quantum level. 

The paper is structured as follows: In Section [IT] we 
briefly review some concepts of QET, introduce the sym- 
metric logarithmic derivative and illustrate the quantum 
Cramer-Rao bound. We also review the notion of dis- 
tance for the quantum Ising model. In Section IIIII we 
derive the ultimate quantum limits to the precision of 
coupling constant estimation at zero temperature, both 
for the case of few spins and then in the thermodynamical 
limit. In Section HVl we analyze the effects of tempera- 
ture and derive the scaling properties of QFI. In Section 
IVl we address the measurement of total magnetization 
as estimator of the Hamiltonian parameter and show its 
optimality. Section IVII closes the paper with some con- 
cluding remarks. 



II. PRELIMINARIES 

In this section we recall the basic concepts of QET and 
the metric approach to quantum criticality. specializing 
them to the one-dimensional Ising model in transverse 
field. 



A. Quantum Estimation Theory 

An estimation problem consists in inferring the value of 
a parameter A by measuring a related quantity X. The 
solution of the problem amounts to find an estimator 
A = A(a;i,a;2, . . .), i.e. a real function of the measure- 
ments outcomes {xk} to the parameters space. Classi- 
cally, the variance Var(A) = -^[A^] — E[X]'^ of any unbi- 
ased estimator satisfies the Cramer-Rao theorem 



Var(A) > 



1 



MF{X) ' 



which establishes a lower bound on variance in terms 
of the number of independent measurements M and the 

Fisher Information (FI) F{X) = E [(^a logp(a;|A)) 



F{\) = Y^p{x\\)[dx\ogp{x\\)Y 



(1) 



p(x|A) being the conditional probability of obtaining the 
value X when the parameter has the value A. When quan- 
tum systems arc involved p(a;|A) — Tv[q\Px\^ {Px^ be- 
ing the probability operator-valued measure (POVM) de- 
scribing the measurement. A quantum estimation prob- 
lem thus corresponds to a quantum statistical model, i.e. 
a set of quantum states p\ labeled by the parameter of in- 
terest, with the mapping A — > pA providing a coordinate 
system. Upon introducing the Symmetric Logarithmic 



Derivative (SLD) Aa as the set of operators satisfying 
the equation 



9a Pa 



1 r 



Aapa + paAa 



we can rewrite the FI as 



FW - E 



Ke{Tr[pxKxPx]f 
Tr [paPx] 



(2) 



(3) 



Then one can prove [ij, ^M that F{\) is upper bounded 
by the Quantum Fisher Information 



F[\) < G(A) EE Tr [paA^] 



(4) 



In turn, the ultimate limit to precision is given by the 
quantum Cramer-Rao theorem (QCR) 



Var(A) > 



1 



MG(A) ' 



which provides a measurement-independent lower bound 
for the variance which is attainable upon measuring a 
POVM built with the cigcnprojectors of the SLD. 

In the following we will consider the quantum statis- 
tical model defined by the set of Gibbs thermal states 
PA = Z-^e-'^"'^^\ {Z = Tr[e-^^(^)]) associated with a 
family H{X) of many-body Hamiltonians where A is the 
coupling constant we wish to estimate. The relevant ob- 
servation at this point is that the Bures distance d^(p, p') 
between quantum states at nearby points in parameter 
space may be written as 

dsl = 4(pA,PA+dA) - gxdX" = \GiX)dX' (5) 

where G(A) is the QFI defined in Eq. ^ and g\ is Bures 
metric tensor given by [22| 



5A 



IE 



\{Tl'm\dxgx\''Pn)\' 
Pn + Pm 



(6) 



\^n) being the eigenvectors of pA = J2nPn\i'n){i'n\. In 
other words, maxima of the Bures metric, e.g. the di- 
vergence occurring at QPTs [8|, correspond to optimal 
estimation working regimes. In the following wc will sys- 
tematically seek for maxima of Bures metric (QFI). In the 
thermodynamic limit those occur at the critical points 
3 , whereas at finite size the maxima of the QFI define 
pseudo-critical points which scale to the actual critical 
points as the size goes to infinity. 



B. Quantum Ising model 

Wc are interested in systems which undergo a zero- 
temperature quantum phase transition and consider a 
paradigmatic example, the one-dimensional quantum 



Ising model of size L with transverse field. The model is 
defined by the Hamiltonian 



h = -jY^ 



f^fcCTfc+l 



k=l 



fc=i 



(7) 



where the af , are Pauli operators and we assume periodic 
boundary conditions cr£_|_x = erf unless stated otherwise. 
As the temperature and the field h are varied one may 
identify different physical regions. At zero temperature, 
the system undergoes a QPT for h = J. For h < J 
the system is in an ordered phase whereas for h > J 
the field dominates, and the system is in a paramagnetic 
state. For temperature T ^ A, A = \J — h\ the system 
behaves quasi-classically, whereas for T 3> A quantum 
effects dominate. The Hamiltonian ([7]) can be exactly 
diagonalized by a Bogoliubov transformation, leading to 



^-Ea^ 



vivk 



1 



(8) 



fc>0 



where A^ denotes the one particle energies and rjk the 
fermion annihilation operator, A^ = -y/e^ + A^, Afc == 
Jsin(A:), €k = (Jcos(fc) + h). Strictly speaking, Eq. ^ 
holds in the sector with even number of fermions. In 
this case, periodic boundary conditions on the spins in- 
duce antiperiodic BC's on the fermions and the momenta 
satisfy k = ^ "^ ''^ . In the sector with odd number of 
particles, instead, one has k ~ LJUI. and one must care- 
fully treat excitations ai k — and k = tt. In any case, 
the ground state of ([7]) belongs to the even sector so that, 
at zero temperature we can use Eq. (|8]) for any finite L. 
At positive temperature we will be primarily interested in 
large system sizes and therefore we can neglect boundary 
terms in the Hamiltonian and use Eq. ([8]) in the whole 
Fock space. For small L we will diagonalize explicitly the 
Hamiltonian ([7]), without resorting to ([5]). 

The QFI for the parameter J may be evaluated starting 
from Eq. ([6|) arriving at 



G.-E 



(djPnf 



2E li-^nldj^. 



,2 jPn -Pmf 
Pn + Pm 



(9) 

which, given Ej ~ X^fc^-fcAfe, where the n^'s are the 
fermion occupation numbers, may be written as [9| 



GjiJ,h,p) = ^E 



4 ^ cosh' (/3Afe/2) 

Ecosh (/3Ai.) — 1 ,„ ,9 / N 
—Zt77^t^idj^>^)- (10) 



cosh(/3Afe) 



where 1?^ — arctan-^. Since the QFI is proportional to 
the Bures metric one may exploit the results derived for 
the Bures metric, which we recall here: 



• At zero temperature, in the off-critical region (the 
thermodynamic limit) L 3> ^, where ^ is the sys- 
tem correlation length, the Bures metric behaves 



as g\ '-^ L\\ — \c\~ " close to the critical point. 
Here L is the size of the system and A^ is related 
to the critical exponents of the transition, which 
turns out to be one for our system [3, Q ■ 

• In the quasi-critical region ^^ L, g\ scales as gx ^ 
L" with a = 1 -I- Ag/i/ where i^ is the correlation 
length critical exponent i.e. ^ ^ |A — Ad^"^. It 
turns out that a = 2 quasi-free fermionic models, 
as the one considered here. 

• At regular points {i.e. not-critical) the Bures met- 
ric is extensive i.e. g\^ L. 

• As the temperature is turned on, as long as it is 
small but larger than the energy gap of the sys- 
tem, quantum-critical effects dominate. In this 
region T ^ A, thermodynamic quantities scales 
algebraically with the temperature and one has 
Qx ~ T"'^, with /3 > 0. For the Ising model /3 = 1 



In the following, we will exploit the dramatic increase 
of the QFI that one experiences in the critical regions, 
to improve the ultimate quantum limit achievable in the 
estimation of the coupling parameter. 



III. QET AT ZERO TEMPERATURE 

In this section we begin to test the idea of estimating 
the coupling constant J of the Ising model by finding the 
maximum of QFI at zero temperature, where the system 
is in the ground state. At first we consider systems made 
of few spins and then we address the thermodynamic 
limit. 



A. Small L 

We start with the case of L = 2, 3 and 4 in Eq.®. The 
QFI is obtained from Eq. ^ by explicit diagonalization 
of the Ising Hamiltonian where p„ — e~'^^"/.Z, En and 
IV'n) are the eigenvalues and eigenvectors of H. For ex- 
ample, for i = 2 we have _E„ — ±2 J, ±2v'^' + h'^ and 
Z ^ 2cosh(2/3J) -h 2cosh(2/3\/j2T7?). For T = one 
gets 



G,/(J,/i,0) 



Gj(J,/i,0) 



h- 


2 

T '^ 


{h' + 


J2)2' ^ - 
3/l2 


4(/l2- 


+ 4/^2 J2 + J4) 



L = 3 



GAJ,h,0) = ^ ^•^(^4 + ;4)2' " ^ L = A. (11) 

Maxima of Gj are obtained for h* — J for L — 2,3,4. 
Actually, this is true for any L, see also the next Sec- 
tion, and the pseudo-critical point h* which maximizes 
Gj^ turns out to be independent of L and equal to the 



true critical point, he = J, VL. At its maximum Gj 
goes like 1/J^ and the ultimate lower bound to precision 
(variance) of any quantum estimator of J scales as J^ . 



We observe that supcrcxtensivc behavior of the QFI in 
the quasi-critical region around the QPT, Gj ~ L^, im- 
plies that the estimation accuracy scales like L~^ at the 
critical points, while it goes like L~^ at regular points. 



B. Large L 

In the following we discuss the QFI for a system of 
size L. We analyze the behavior of Gj near the critical 
region at T = 0. Taking the hmit T ^ in Eg. pH)) . 
the classical elements of the Burcs metric, which depends 
only on thermal fluctuations, vanishes due to the factor of 
(cosh(/3Afc/2))~^. Therefore, at zero temperature, only 
the nonclassical part of Eq. ([TO)) survives and one obtains 



(12) 



Gj = Y.^dj^k? 

k 



where djdk = ,+(^^\,,,y. {dj^) = ^^. Since we 
are in the ground state, the allowed quasi-momenta are 



(2Tt+l)Tr 



with 



0,. 



Gj = 



E 

k 



L/2 — 1. Explicitly we have 
/i2 sin(fc)2 



K 



(13) 



We are interested in the behavior of the QFI in the quasi- 
critical region £, ^ L. In the Ising model i/ = 1 so the 
critical region is described by small values of the scaling 
variable z = L(/i — J) ~ i/C; that is 2; « 0. Conversely 
the off-critical region is given by z — > 00. We substitute 
h — J + z/L in Ea. p3p and expand around z = to 
obtain the scaling of Gj in the quasi-critical regime 

G -V (J+f)^sin^(fc„) ^^ 

' r [|^+4J(J+f)sin2(fc„/2)P Z.^'="^^ 

(14) 
Since dzf{0) = 0, the maximum of Gj is always at z = 
for all values of L, in turn, the pseudo-critical point is 
h*^ = J = he Vi. As already noticed previously, the 
statement /i^ — he is peculiar to this particular situa- 
tion. For instance, introducing an anisotropy 7 so as to 
turn the Ising model into the anisotropic XY model, the 
pseudo-critical point gets shifted and one recovers the 



general situation h\ 
order one obtains 



^4J2 



0{L- 



cot"(fc„/2) 1- 



Going to second 



2J2L2 



sin2(fc„/2) 



)+0{z'). (15) 



Using Euler-Maclaurin formula [10[ we get 



^HsJ^ -3837^ -875 + °<^°>^ («) 



This shows explicitly that aX h ~ J the Fisher informa- 
tion has a maximum and there it behaves as 

8J2 



Gj{L,T^O,h* ^J) 



0{L). 



(17) 



C. Signal-to-noise ratio 

Notice that, in assessing the estimability of a param- 
eter A, the quantity to be considered is the quantum 
signal-to-noise ratio (QSNR) given by Q{\) = \^G{X) 
which takes into account of the scaling of the variance 
and the mean value of a parameter rather than its ab- 
solute value. We say that a parameter A is effectively 
estimable when the corresponding Q(A) is large and that 
to a diverging QFI corresponds the optimal estimability. 
In both cases of few and many spins, at the critical point 
the QFI goes like 1/ J^, this means that Q{J) is indepen- 
dent on J and one can estimate large as well as small 
values of parameters without loss of precision. 



IV. QET AT FINITE TEMPERATURE 

In this section we consider the estimation the coupling 
constant J at finite temperature. We first discuss in some 
detail the short chains with L = 2,3,4 and then we treat 
the case L ^ 1. 



A. Small L 

As a warm-up let us first focus on the simplest, L = 2 
case. A first step in the computation of the SLD for 
two qubit is to find the SLD in the single qubit case. 
Consider a system with " Hamiltonian" H — a- cr in the 
state p = e~^'^Z~^ where Z = Tre""**^, and the three- 
component vector a depends on parameter J. The SLD 
relative to this state turns out to be 



A = — tanh (a) idjk ■ cr) 



1 + tanh (a) - 2 tanh (a) (dja) (a • cr) . (18) 



where a is the modulus of a and a = a/a. Now 
note that the Hamiltonian Q for L = 2 (with PBC), 
has the following block-diagonal form in the basis 



H^-2f3 



Ja'^ + ha"" 
Jcr-^ 



(19) 



We can then apply formula ([TB]) in each subspace to ob- 
tain the full SLD. After some algebra one realizes that 
the SLD has the following form 



A = cicr^ 



' -f C2Cr« ® cr^ -f- C3 (cr^ (g) I + I ® cr^) , (20) 



where ci^2,3 are constants which depend on /3, J, and h. 
When the temperature is sent to zero the above expres- 
sion becomes 



A 



T=0 — 



/iK 



2(j2+/l2)3/2 



(21) 



We see that, aheady in the simple two-qubit case, the 
SLD is a comphcated operator both at positive and at 
zero temperature. More involved expressions are ob- 
tained L = 3,4 and larger. We do not report here the 
analytic expression of the corresponding QFIs Gj for 
L = 2,3,4 since they are a bit involved. Rather, in or- 
der to assess estimation precision at finite temperature 
and compare it to that at T = 0, we consider the ratio 
-fj = Gj{f], J, h)/Gj{oo, J, h), for some fixed values of J 
and illustrate its behavior by means of few plots. 



This is illustrated for L = 2 in Fig. [2l where we re- 
port the log-linear plot of Gj{f3) as a function of h for 
different values of J and /3. For high temperature the 
maxima are localte at a field value close to zero, whereas 
for decreasing temperature they move towards h* = J. 



Log Gj{y8) 




Figure 2: Logarithmic plot of Gj{P) vs h for (from top to 
bottom) J = 0.1, 0.5, 1, 2,5, 8 and /3 = 1 (left), /? = 10 
(right). 




Figure 1: The ratio 7j as a function of h for L = 2, 3, 4 (from 
top to bottom) and J = 0.5 (on the left), J = 5 (on the 
right). The colors refer to different values of /3 = 1 (red), 
/3 = 10 (green), (3 = 100 (purple), /3 = 1000 (blue). 

As it is apparent from Fig. [T] for small h the ratio is 
smaller than 1, i.e estimation of J is more precise at zero 
temperature, whereas for increasing h a finite tempera- 
ture may be preferable. In turn, for any value of J and /?, 
there is a field value that makes finite temperature con- 
venient: this is true also for low temperature as proved 
by the presence of a global maximum for small h, besides 
the local maximum at h = J. For /3 — > oo the maxima at 
small h disappears and we recover the zero temperature 
results. Notice that, in view of Eqs. pTj) . the ratio 7j 
is proportional to the QSNR. Besides, since the maxima 
of 7j vary with /3, we conclude that the optimal field 
h* which maximizes Gj{(3) varies with the temperature. 



B. Large L 

At positive temperature and L large, the sums in equa- 
tion (jlOp may be replaced by L J dk. The quantity 
Gj = Gj/L is always convergent, and the convergence 
rate is exponentialy fast in L in the (renormalized clas- 
sical) region T ^ A whereas is effectively only algebraic 
when T 3> A (the quantum-critical region). Thus, up 
to contribution vanishing with L, Gj = G]j + G^ is a 
bounded function of its arguments as long as T > 0, 
given by 



G\ = ^ 



(3' 



dk 



{J + h cos {k)Y 



ttJo cosh' (/3Afe/2) 



A? 



Gi = 



1 

2^ 



dk 



cosh(/3Afc)- lh'^sui{ky 
cosh(/3Afc) A4 



(22) 



(23) 



For any T > the function Gj has a cusp in h = J, 
where it achieves its maximum value. Changing variable 
from momentum to energy, the integrals above can be 
approximately evaluated in the quantum critical region 
/3 I J — /i| ^ 1 (actually we also require low temperature, 
i.e. /3 I J + /i| > 1). The result is 



G'j = 



G'j 



% 


(3) T 


8 
C 


K .P\J+h\ 
\J + h\ 1 



TJ2 



8J2 



O (T°) 

f 0(T) 



(24) 
(25) 



where C is Catalan's constant C — 0.915 and the Riemann 
Zeta-function gives C(3) = 1.202. 

In summary, for large sizes and at positive temper- 
ature, the maximum of the QFI as a function of h is 
always located ai h = J for all values of J, T. At the 
maximum, the QFI is approximately given by 



Gj 



2C L 



(26) 



As a consequence, the QSNR scales as Qj ~ JL/T, 
in other words, at finite temperature, the estimation of 
small values of the coupling constant is unavoidably less 
precise than the estimation of large values. As expected, 
large L and/or low temperature improve the precision of 
estimation. 



V. PRACTICAL IMPLEMENTATIONS 

The SLD represents an optimal measurement, i.e. the 
corresponding Fisher information is equal to the QFI. 
However, as we have seen [see e.g. Eq. ([2T|)]. generally the 
SLD does not correspond to an observable whose mea- 
surement can be easily implemented in practice. There- 
fore, in this section, we consider the total magnetiza- 
tion Mz ~ -^ X^i '^i^ ^^ ^ feasible and natural measure- 
ment to be performed on the system in order to esti- 
mate the coupling J. We assume that the system is at 
thermal equilibrium, p ~ Z^^e^^^ , and consider short 
chains L = 2,3,4. We illustrate the procedure in de- 
tail for the simplest L = 2 case. Upon measuring Mz, 
the possible outcomes are m = {1,0,— 1} with eigenpro- 
jcctors P™ given by Pi = |00)(00|, P_i = |11)(11|, and 
Po = |10)(10| + |01)(01|. The corresponding probabilities 
p(m\J) = Tr{pPjn) are given by 



p(±l| J) = 



3h(2/3VJ2+7?) 



p{0\J) 



2 [cosh(2/3 J) + cosh(2/3\/J2 + /i2)] 

(l ± h{J^ + h^)-^/^ tanh(2/3v/ J2 + h^) 

cosh(2/3J) 
'cosh(2/3J) + cosh(2/3\/j2T7?) 



X (27) 



(28) 



The FI is then obtained by substituting p(m\J) into 
Eq. ([T]). The resulting expression provides a bound for 
the variance of any estimator of J based on M measure- 
ments of magnetization: Var( J) > l/MFj. 

Fj(r)/Gj(h*) 
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Figure 3: The ratio Fj{P, J, h~)/Gj{P, J, h*) as a function of 
J for L = 2 (red line), i = 3 (green line), and L = 4 (blue 
line). The bottom group of lines is for /3 = 3, whereas the top 
group is for /3 = 10. 

The Braunstein-Caves inequality says that the FI of 
any measurement Fj is upper bounded by the QFI Gj. 



For the magnetization this is illustrated in Fig. [3l 
where we plot of the ratio Fj{(3, J,h~)/Gj{P,J,h*) for 
L = 2, 3, 4, h~ being the field maximizing the FI. Notice 
that for increasing J the FI of the magnetization saturate 
to the QFI, i.e. magnetization measurement becomes op- 
timal. The saturation is faster for lower temperature (we 
report the ratio for /? = 3 and (3 = 10). Notice also that 
for low temperature the dependence of the ratio on the 
size L almost disappears. In summary: for any temper- 
ature there is a threshold value for J, above which the 
measurement of the magnetization is optimal for the es- 
timation of J itself. This threshold value decreases with 
temperature, and for zero temperature magnetization is 
optimal for any J. Indeed, after explicit calculation of 
Eq. dl]) for i = 2, 3, 4 we found that, in the limit T -^ 0, 
Fj{h,T = 0) = Gj{h,T ^0),i.e. the FI of the magneti- 
zation is equal to the QFI. In other words the estimation 
based on magnetization may achieve the ultimate bound 
to precision imposed by quantum mechanics. Besides, 
at finite temperature, despite the fact that the equal- 
ity does not hold exactly, Fj is only slightly grater than 
Gj almost in the whole parameter range {J,T). This 
may be also seen in the behavior of Fj versus temper- 
ature: the ratio Sj = Fj{f3,J,h)/Fj{oo,J,h) at fixed J 
may be greater than 1 for some values of the magnetic 
field, namely magnetization measurements may be more 
precise at finite T, as it happens for the optimal mea- 
surement with precision bounded by the QFI. Of course, 
for T ^ 0, Sj -^ 1. 

Overall, we conclude that the magnetization Af^ is a 
good candidate for nearly optimal estimation. Of course 
we still need an efficient estimator, that is an estima- 
tor actually saturating the (classical) Cramer- Rao bound. 
To this aim we employ a Bayesian analysis, since Bayes 
estimators are known to be asymptotically efficient [23|, 
i.e. VaT{JB) = l/MFj, for M > 1, Js being the 
Bayesian estimator (see below). According to the Bayes 
rule, given a set of outcomes {m} from M independent 
measurements of the magnetization, the a-posteriori dis- 
tribution for the parameter J is given by p{J\{m}) = 
l/iV]^j^jP(r7i| J)"™ where A'^ is a normalization constant 
and rim is the number of measurements with outcome 
m. Bayes estimator is the mean Jb = J dJ J p{J\{m}) 
of the a-posteriori distribution and precision is quanti- 
fied by the corresponding variance. In the asymptotic 
limit of many measurements M ^ 1, n^ — > Mp{m\J*) 
where J* is the true value of the parameter to be es- 
timated and the a posteriori distribution rewrites as 
Pa{J\{m}) = l/A^^„^exp[Mp(m|J*)logp(m|J)]. In or- 
der to check the actual meaning of " asymptotic" wc have 
performed a set of Monte Carlo simulated experiments of 
the whole measurement process. In Fig. |4l wc report the 
result of Monte Carlo simulated experiments of magne- 
tization measurements for ,7 = 3 and (3 = 1. The black 
dots represent the mean variance of the Bayes estima- 
tor Jb averaged on 20 sets each of 500 measurements. 
The blue line is the corresponding variance evaluated us- 
ing the asymptotic a-posteriori distribution, whereas the 
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Figure 4; Estimation of the coupling constant by magneti- 
zation measurement and Bayesian analysis. The plot is for 
L — 2 and P — 1 with true value of the coupling equal to 
J = 3. We report the variance of the Bayesian estimator 
for Monte Carlo simulated experiments (M=500, black dots), 
the corresponding variance evaluated using the asymptotic a- 
posteriori distribution (blue line), and the (classical) Cramer- 
Rao bound (green line). 



green line is the Cramer-Rao bound {MFj)^^ . The plot 
shows that Bayes estimator is indeed asymptotically ef- 
ficient and that already with a few hundreds of measure- 
ments one may achieve the ultimate precision. Overall, 
putting this result together with the fact Fj ~ Gj (see 
Fig. [3]) we conclude that the measurement of the total 
magnetization of the system provides a nearly optimal 
and feasible measurement (at any (3) to estimate the cou- 
pling of the one-dimensional quantum Ising model. 



VI. CONCLUSIONS 

The coupling constant of a many-body Hamiltonian 
is not an observable quantity and we have to solve a 
quantum statistical model to evaluate the bounds to its 
estimation precision. This has fundamental implications 
since it corresponds to find the ultimate limits imposed 
by quantum mechanics to the distinguishability of dif- 
ferent states of matter. In this paper we exploited the 
equivalence between the quantum Fisher metric and the 
(ground or thermal) Bures metric and all the results re- 
cently obtained for the latter. Specifically at zero tem- 
perature, the Bures metric scales with the system size 
L at regular points whereas it can increases as L^ at or 
in the vicinity of quantum critical point. A similar en- 
hancement takes place when temperature is considered. 
In turn it is possible to exploit this enhancement to dra- 



matically improve the bounds to precision in a quantum 
estimation problem. Let us imagine that an experimenter 
would like to infer the value of a coupling constant of a 
physical system over which he has little or no control. 
Reasonably the experimenter has good control over the 
external fields he can apply to the system. The idea is 
then to tune the external field to a value close to the 
quantum critical point. At this value of the couplings, 
an improvement of order of L can be achieved in the pre- 
cision of the estimation of the unknown coupling. To 
test these ideas in practice, we have worked out in detail 
a specific example, the ID quantum Ising model. This 
model provides us with all the ingredients we need, a 
coupling constant J, an external field ft., and a quantum 
critical point at h ~ J. The main accomplishments of 
our analysis are: i) At zero temperature we evaluated 
the precision in the estimation of the coupling, exactly 
for short chains oi L ~ 2,3,4 sites and asymptotically for 
large L. We found that the optimal estimation is possible 
at values of the field exactly equal to the critical point, 
independently of L. For large L we indeed observe a 1/L 
enhancement of precision, and a quantum signal-to-noise 
ratio independent of the coupling, ii) At positive tem- 
perature the optimal value of the field is again given by 
the critical value when the system size is large or the 
temperature is low. In the other working regimes the op- 
timal field maximizing the quantum Fisher information, 
defines a set of pseudo-critical points, and the optimal 
precision scales as TJ/L. iii) We obtained the optimal 
observable for estimation in terms of the symmetric log- 
arithmic derivative and showed that already in the case 
L = 2 it does not correspond to an easily implementable 
measurement, iv) We have shown that for small L the 
measurement of the total magnetization allows to achieve 
ultimate precision. Using Monte Carlo simulated exper- 
iments and Bayesian analysis we proved that this is pos- 
sible already after a limited number of measurements of 
the order of few hundreds. We conjecture that this may 
be true for any L; work along this line in progress. 

Overall, we found that criticality is a resource for 
precise characterization of interacting quantum systems 
(e.g. a quantum register), and may represent a relevant 
tool for the development of integrated quantum networks. 
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